Reduction of coating thermal noise by using an etalon 
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Reduction of coating thermal noise is a key issue in precise measurements with an optical inter- 
ferometer. A good example of such a measurement device is a gravitational-wave detector, where 
each mirror is coated by a few tens of quarter-wavelength dielectric layers to achieve high reflectivity 
while the thermal-noise level increases with the number of layers. One way to realize the reduction 
of coating thermal noise, recently proposed by Khalili, is the mechanical separation of the first few 
layers from the rest so that a major part of the fluctuations contributes only little to the phase 
shift of the reflected light. Using an etalon, a Fabry-Perot optical resonator of a monolithic cavity, 
with a few coating layers on the front and signiflcantly more on the back surface is a way to realize 
such a system without too much complexity, and in this paper we perform a thermal-noise analysis 
of an etalon using the Fluctuation-dissipation theorem with probes on both sides of a finite-size 
cylindrical mirror. 
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1. OVERVIEW 

Brownian fluctuation of multi-layer coatings on a 
mirror is a dominant noise source in interferometric 
gravitational-wave detectors. One way to reduce coat- 
ing thermal noise is mechanical separation of the first 
few layers from the rest of the coatings [H0. This so- 
called Khalili-cavity will be useful in precision measure- 
ments like the sub-SQL measurement experiment Q or 
a future-generation gravitational-wave detector like Ein- 
stein Telescope Replacing a single mirror by a 
cavity will, however, increase the complexity of the sys- 
tem, and a way to ease this problem is to use an etalon 
instead of the double-mirror cavity. The front surface 
and the back surface of the etalon are not mechanically 
separated, but as we will show in this paper, the separa- 
tion is good enough to reduce the noise level of the whole 
system better than a single conventional mirror. 

Figure [1] shows our model. The test mass is a cylin- 
drical mirror made of silica and its motion is probed by 
an axisymmetrie Gaussian beam. Both sides of the mir- 
ror are coated by silica-tantala doublets and the distance 
of the two surfaces is controlled to be anti-resonant for 
the carrier light utilizing the temperature dependence of 
the refractive index of the substrate. A larger fraction 
of the light is reflected back by the front surface but the 
thermal-noise level on the front surface is lower than that 
of the back surface because of the fewer coating layers. 
Each surface of the etalon is probed less than a fully re- 
flective surface of a conventional mirror; the front surface 
of the etalon is probed by |ei| < 1, and the back surface 
is probed by |e2| < 1. 

Total coating thermal noise of the etalon is the sum 
of noise on the front surface and noise on the back sur- 
face, which are caused by the thermal energy stored in 
the coating layers. In fact, it is not only the energy 
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FIG. f : Khalili etalon with a few coating layers on the front 
surface and more layers on the back surface. The light is 
anti-resonant in the etalon so that a larger fraction of the 
light incident onto the etalon is reflected by the front surface. 
Mechanical losses of the back-surface coatings cause phase 
noise on the probe beam not only via thermal motion of the 
back surface, which is measured less, but also via thermal 
motion of the front surface as the two surfaces are connected 
by the substrate with flnite thickness. 



stored in the front (back) surface coatings that causes 
the phase fluctuation on the beam reflected at the front 
(back) surface of the etalon. The thermal energy stored 
in the coatings of the other surface contributes to the 
phase fluctuation on the beam via the substrate. For the 
etalon, we should thus consider four elements to calcu- 
late total coating thermal noise. One is the fluctuation 
of the front surface caused by the thermal energy in the 
front surface coatings. This process is the same as that of 
coating thermal noise of a conventional mirror. Second is 
the fluctuation of the front surface caused by the thermal 
energy in the back surface coatings that is mechanically 
transferred to the front surface. The other two are the 
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fluctuation of the back surface coatings caused by the 
thermal energy in the front surface coatings and that in 
the back surface coatings, which are probed by the beam 
transmitting through the substrate or in other words op- 
tically transferred to the front surface. The fluctuations 
caused by the thermal energy in the coatings of a same 
surface are correlated. 

According to the Fluctuation-dissipation theorem, the 
fluctuation caused by the thermal energy is given by cal- 
culating dissipation of the elastic energy caused by an 
imaginary force that is intentionally applied to the mea- 
sured surface the case of a conventional mirror, 
a single imaginary force is applied to the reflective front 
surface, and the elastic energy is integrated over the coat- 
ings on the surface. In the case of an etalon, an imagi- 
nary force shall be applied to each reflective plane with 
a weight coefScient ej, which is given from the response 
from the motion of each surface to the phase change of 
the reflected light. The elastic energy is then integrated 
over the coatings on both surfaces. The thermal-noise 
level is proportional to the dissipated power, which is 
given by the elastic energy multiplied by the mechanical 
loss angle. 

To calculate the clastic energy of the etalon, we should 
solve the elastic equation of the finite-size cylinder. As is 
done in Ref. , we use Levin's method with one probe 
beam on the front and another probe beam on the back 
surface. The way to extend Levin's method to a finite 
size mirror has been shown by Bondu et al. [lo| 11 



and we follow the same way with slightly different bound- 
ary conditions using the two probe beams; a similar ap- 
proach has been introduced without the derivation in 
Ref. [m to evaluate thermal noise of a monolithic cav- 
ity. We will show the semi-analytical derivation for the 
model with the two probe beams in Sec. [51 In fact, for 
this model, the front surface coatings are regarded to be 
so thin that the difference between the front side and the 
rear side of the coating layers is negligible - thin-layer 
model. A closer look into the optical behavior reveals that 
the light transmitting through the front surface coatings 
circulates inside the etalon with reflecting on the back 
surface coatings and on the rear side of the front sur- 
face coatings. In Sec. [3l we take the difference of the two 
sides of the front surface coatings into account - thick- 
layer model. Section|4]shows the result of the calculation. 
We have verified the result by numerical calculation using 
a finite-element analysis code, which is shown in Sec. [5] 
with other discussions. 
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summed (see Sec. 15.41 for details). The reflectivity of the 
approximated mono-layer coating is that of the multi- 
layer coatings. 

In the thin-layer model, we assume that the light is 
always reflected at the front side of the coatings. This is 
the case for a conventional mirror where the light trans- 
mitting through the coatings never comes back. This is 
also true for the back surface coatings. The question is 
about the front surface coatings, but let us ignore the 
difference for the simplicity. This simplification is valid 
as the anti-resonant etalon accommodates little light cir- 
culating in it. 



2.1. Probe force amplitudes ej 

We can use the quasi-static approximation for a short 
cavity, which means that displacements of mirrors are 
sufficiently slow compared with the relaxation rate of 
the cavity. The optical path length between the surfaces 
can be assumed to be an odd multiple of the quarter- 
wavelength in the zeroth order. Phase of the light re- 
flected by the cavity as a compound mirror changes ac- 
cording to the differential motion of the two mirrors 
dx2f — Sxif. In addition, the phase shift due to the mo- 
tion of the front mirror from its initial position should 
be taken into account. Note that the former motion is 
probed inside the substrate, thus the phase shift is pro- 
portional to the refractive index of the silica substrate 
Ug. The reflected light field Eout can be given with the 
input field Ein as 
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E,n l + rir2e2*« 

9 = kons{Sx2f - Sxif) , 



(1) 



rj is the amplitude reflectivity of the front/back surface, 
and fco is the wave number. Expanding Eq. ([T|) into a 
series over Sxn and Sx2f, and keeping linear terms, we 
obtain 
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Ej. 



ri + r2 
1 + rir2 



2iko {eiSxii -I- €26x21) 
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with 



ei 



£2 



ri+r2 nsr2{l-rl) 



l-\-rir2 (l + rir2)2 
nsr2il - rl) 
(l + rir2)2 



(3) 



These ej are used as weight coefficients of the imaginary 
forces applied to the surfaces of the etalon. 
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Figure [T] shows our model. As is done in Ref. 
the multi-layer silica-tantala doublets are approximated 
as a single thick layer of silica or tantala. The dissipa- 
tions arc calculated for silica and tantala to be square- 



2.2. Elastic motion of the substrate 

The Fluctuation-dissipation theorem tells us that the 
power spectrum of thermal noise at the radial frequency 
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n is given by 



TTA 
~U(I) , 



(4) 



where is the Boltzmann constant, T is the tempera- 
ture, U is the maximum elastic energy that can be gen- 
erated by the imaginary force Fq, and is the loss angle. 

While the conventional method applies a single imagi- 
nary force on the only reflective surface of a mirror, our 
new method shall apply imaginary forces on both sides 
of the etalon. The weighting coefficients of the imaginary 
force are ei and £2 that we have derived above. 

The elastic energy is given by the product of the strain 
tensor Eij and the stress tensor Tij, integrated over the 
volume of interest: 
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We use a cylindrical coordinate system along the z axis; r 
is the distance from the z axis and ip is the angle around 
the z axis. Our mass is a cylinder with radius of a and 
thickness of h. The strain tensor elements of the cylin- 
der with the axisymmetric pressure are expressed by the 
displacement vectors Ur and Uz as follows: 



Er 



dur 
or 



r 



dz 



Erz — — 



dUr 



dr 



(6) 



and the stress tensor elements are as follows: 

T„. = {\ + 2n)Err + \{E^'4,+Ezz) 

T^tp = (A + 2fl)E^^ + X{Ezz + Err) 

Tzz = {X + 2fl)Ezz + HErr + E^,^) 

Trz — 2fJ,Erz . 

Here A and fi are the Lame coefficients: 
Yu Y 



A = 



(l + z/)(l-2i/) ' 



2{1 + iy) ' 



(7) 
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with Y the Young^ modulus and v the Poisson's ratio. 
According to Ref . [9[ , displacement vectors of a finite-size 
cylinder are described with the Bessel Functions: 

Ur{r,z) = ^^Am{z)Ji{kr) + SUr , (9) 

m 

Uz{r,z) = Y,Bm{z)Jo{kr) + Suz . (10) 



Am and _B,„ are in the same form as in Ref. 

Am{z) = i,ne~''^ + 5me^' 



kz X + ^ 
X + 2ii 

X + 3fi 



2(A + 2/i) 



-kz 



-kz 



kz-\ 



(11) 



X + 3fi 

2{X + 2n) 

kz X + u 

a 

2 A + 2/.t 



5ur and 5uz are given as 



5ur 
Suz 



air + a2rz , 

+ 62r^ + b^z 



(12) 

(13) 
(14) 



The coefficients 7™, 5„i, ai, 02, bi, 62, and 63 

are determined with the following boundary conditions: 



Tzzir,0) ^ -eiFop{r) 
Tzz{r,h) = e2Fop{r) , 
Trz{r,0)^0, 
Trz{r,h)=0 , 
Trz{a,z) ^ , 

Trr{a, z) — . 



(i) 

(ii) 

(ill) 

(iv) 

(v) 

(vi) (15) 



We have ej on the right side of the boundary conditions 
(i), (ii); the thermal- noise level of a single mirror will be 
given if we put ei = 1 and £2 ~ 0. The beam profile p{r) 
is expanded with Bessel functions: 



P(^) = y^^PrnJajkmr) +po , 

rn 

exp [-kf^w'^/8] 



■na^Jl[k„ia) 



(16) 

(17) 
(18) 



Here km k) arc the zeros of the first-order Bessel func- 
tion, divided by the mirror radius a; namely Ji(fca) = 0. 
The Pn term is missed in Ref. Q and is corrected in 
Ref. [lOl . Boundary condition (i) reads 



f k^i[{am - Pm) + 2(7m -I- 5m)] = ^iFaPr, 
\ (A + 2/i)63 + 2Aai = -eiF^po . 

Boundary condition (ii) reads 



(19) 



fc/i 



A + ^j 



I3me'"') + 2{^me'^'' +5me^'') 



A + 2^i 

^ (A + 2n){2bih + 63) + 2A(ai + 02/1) = e2^oPo 



(20) 
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Boundary condition (iii) reads 



2kfi 



{am + /3m) + 2 (7m - 6m] 



Boundary condition (iv) reads 



. (21) 



Boundary condition (v) reads 

(as + 262)0 = 0. 



(23) 



2fc/i 



kh\ 



X + 2iJ, 



= 



Boundary condition (vi) will be almost satisfied by min- 
imizing / = J^'' T^^{z)dz. 

Trr = e(z) + Co + CiZ , (24) 

e = [(A + 2ii)kAm + XB'J Jo{ka) , (25) 

Co = 2(A + ^)ai + A63 , (26) 

ci = 2(A + Ai)a2 +2A61 . (27) 

(22) From Eqs. (Il9l)-(l22l) we obtain 



FoPmjX + 2^l) ei(l - Q + 2khQ) + £27^(1 + 2fcfe - Q) 
kii{X + ii) (1 - g)2 - 4fc2/i2g 

Fop,„(A + 2/^)Q £1(1 - Q + 2fcfe) + e2lVQ{l + 2fcfeQ - Q) 
k^i{X + fi) (1 - Q)2 - 4fc2/i2Q 

Fop™ £i([2fc2/i2(A + ^) + 2^fc/i]g + ^(1 - Q)) + e^VQiK'^ ^ Q) + kh[{X + /i)(l - Q) + 2/i]) 



2fc/Lt(A + /x) 
FoPmQ ei{[2eh^{X + fi) -2^ikh]- 



(1 - Q)2 - Ak^h^Q 

- Q)) + £2/VQ(-/^(l -Q) + kh[{X + /i)(l - 0) - 2^Q]) 



2fc^(A + ^) 



where Q = exp[— 2A:/i]. Also 61, 62, 63 can be expressed 
by ai and a2 as: 



(1 - g)2 - 4fc2;i2Q 



(28) 



we obtain 



hi 

&2 



X 



1 



-02 



A + 2^ 2(A + 2/i) 
02 
2 ' 

2A 1 



(ei +£2) 



Q{z)dz = , 



(29) 



/ Q{z)zdz = 



£1 + £2 

/c2 



FoPmJo{ka) . 



A + 2/i X + 2n 
Now let us take the derivative of / 



Combining these with Eqs. (|30l) and (|3T]) . we have 



5/ 



1 



— =0^co = --/ e(z)dz--ci, (30) 



E 



6(£l + £2) 



'^0 ^ ^ ' fc2;^2 PoPmJaika) 



Cl 



E 



12(£l + £2) 

k^h^ 



Fopm-h{ka) 



(32) 

(33) 

(34) 
(35) 



(31) 



After some algebra with the boundary conditions (i)-(iv), 

A + 2/i 



and then ai and 02 are given from Eqs. ((26| and ((27l 
Finally, we obtain 



5ur = 



(cor + Cirz) + 
2 



CqZ - 



ClZ 



2^(3A + 2^i) 
X 

~ fi{3X + 2^i) \ '" ' 2 
(A + fJ.)Fopo 



2n{3X + 2pi) 



{ei-iei+e,)^} , 



X + 2fi 2 
-cir 



4/i(3A + 2^) 



Ai(3A + 2/i) 

As is introduced by Harry et al. [3], the boundary 



£lZ- (£1 + £2)77r f + 



AFopo(ei + £2)?-^ 



(36) 



(37) 



2ft. j 4/i(3A + 2/x)/i 
conditions between the substrate and the front-surface 
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coatings are: 



E'„{r 
TLir 



Err{r,0) , 

T„(r,0) = , 
T..(r,0) , 



(38) 



where the elements with a prime (') are for the front- 



surface coatings. Since our probe forces are oriented 
afong the normal of the cylinder's flat faces (i.e. z-axis) 
they induce no shear to the etalon. Consequently, at the 
boundary and in the coatings Erz, E'^^, Trz, and T^^ are 
zero. Since the coatings are thin, we can assume that the 
strain and stress tensor elements are constant in terms of 
z. The strain tensor elements of the front-surface coat- 
ings are then given as 



k r 



E'„ = J2^miO)- Mkr)~.h{k 



X + 2fi 



2/i(3A + 2^)^° ' 2n{3X + 2fi) 



X + 2fi 



2fi{3X + 2nf° ' 2fi{3X + 2fi) ' 



XeiFopo 



e' — 

X' + 2n' 



kAm\0)Jo{kr) + , „ , co ' 



Ai(3A + 2/i) ^i{3X + 2^l) 



(39) 



and the stress tensor elements of the coatings are given 
as 

Trr = (^' + 2/i')-^rr + + -E-L) I 

T^,^ = {X' + 2^')E'^,^, + X'{E'^^+ E[.^) , 

= (A' + 2/i')i?L + ^'{Kr + E'^^) ■ (40) 

Putting these into 

U' = 7Tdi f E'^jT^j^rdT (j = r, V, z) , (41) 
Jo ^ 

and then into Eq. (|4]), we obtain the power spectrum 
of thermal noise originating from the loss in the front- 
surface coatings. 

The boundary conditions between the substrate and 
the back-surface coatings are: 



E':,.{r 
E'Uir 

E'lAr 
T'JAr 



Err{r, h) , 

E.4,^{r,h) , 
ErAr, h)=0 , 
TrAr,h)=0 , 



(42) 



where the elements with a double prime (") are for the 
back-surface coatings. Note that E'^'^, T"^, and r"^ 
zero. Since the light on the back surface is reflected back 
at the boundary of the substrate and the coatings, the 
imaginary force — e2^oJ'('') is applied at the boundary. 
The stress in the z-dircction in the back-surface coatings 
is then zero. This means that the expansion of the back- 
surface coatings does not appear in coating thermal noise. 
The strain tensor elements of the back-surface coatings 
are given as 



i?;; = ^A^(/.)^ 



Jo{kr) - J2{kr 
E;^ = ^ A„(/j)| [jo(fcr) + J2{kr) 



E'' 



-A" 



A" + 2^ 



-^/cyl,„(/i)Jo(fcr) 



A + 2^ 
2/z(3A-h 2^) 

A + 2^ 



A£2-FoPO 

2n{iX + 2ii) ' 
Ae2foPo 

2^(3A + 2/^) ' 2^(3A + 2/^) ' 

A" [(A + 2m)(co + cih) - Aezfopo] 
(A" + 2/i")/i(3A + 2/_i) 



(co + cih) - 
(co + cih) - 



(43) 
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and the stress tensor elements of the eoatings are given 
in the same way as in Eq. (|40)) . Putting these into 



Jo 



r, ^, z) , (44) 



and then into Eq. (jlj, we obtain the power spectrum of 
thermal noise originated from the loss in the back-surface 
coatings. 
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increase much. We apply the imaginary force as usual on 
the rear side of the front surface coatings, namely the 
interface between the coating layers and the substrate. 
It is the stress tensor inside the front surface coatings 
that makes a slight difference from the thin-layer model. 

In the thin-layer model, we applied the imaginary force 
with the weight function ei on the front side of the front 
surface coatings, which imposes thermal noise due to the 
expansion of the coatings in the z-direction by ei. If we 
apply the imaginary force on the rear side of the front sur- 
face coatings, the expansion is not probed by the light, 
which is the case for the back surface coatings. What 
happens in reality in the front surface coatings is inter- 
mediate of the two extreme situations. A fraction of the 
light is reflected on the front side and probes the expan- 
sion, and the rest transmits through the coatings to be 
reflected by the back surface coatings of the etalon and 
nevertheless probes the expansion but in a different way 
due to the difference of the refraction index. 



3.1. Probe force amplitudes tj 
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FIG. 2: Top : Thin-layer model. The light at the front surface 
is always reflected on the front side of the coatings. Bottom : 
Thick-layer model. The light circulating inside the etalon is 
reflected on the rear side of the front surface coatings while 
the incident light is reflected on the front side of the coatings. 
The mathematical formulation with the thin-layer model is 
more simple while the thick-layer model is more realistic. 



3. THICK-LAYER MODEL 

In the thick-layer model, we take into account the fact 
that the light circulating in the etalon probes the rear 
side of the front coating. Figure [2] explains the difference 
between the two models. Imaginary forces are applied at 
three different planes in the thick-layer model since the 
reflection planes are actually different for the incident 
light and the light circulating in the etalon substrate. In 
fact, since there are anyway two boundary conditions on 
Tzz , the complexity to solve the elastic equation does not 



The optical path length between the rear side of the 
front surface and the front side of the back surface is 
kept an odd multiple of the wavelength in the zeroth 
order. The optical thickness of the front surface coatings 
is also kept an odd multiple of the wavelength in the 
zeroth order. Phase of the light reflected by the etalon 
changes according to the motion of each boundary Sxn, 
Sxvo, and 5x2i- The reflected light field Eout is given 
with the input field Ein as 



Ei, 



2ie, 



(1 - rl)r2e 



2ie 



1 



rir2e' 



fconc((5a;ib - faif) , 



(45) 



where ric is the refractive index of the coating layers; the 
index for silica to calculate thermal noise of silica coatings 
and the index for tantala to calculate thermal noise of 
tantala coatings. Expanding Eq. (j45p into a series over 
(5xif , (5xib, and 5x2i, and keeping linear terms, we obtain 



-Bo 



Ei, 



with 



1 + rir2 



- 2iko (ciffait -I- eihSxih + e2f<5a:2f ) , 



(46) 



eif 



eib 



e2f 



ri+r2 nc{l-rl)r2 



1 + rir2 
?ie(l - rl)r2 

1 + 'rir2 
Usjl - rl)r2 
(l + rir2)2 



+ rir2 
ns(l - 



(l + rir2)2 



(47) 
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These ej are combined to be used as weight coeffi- 
cients of the imaginary forces applying to the surfaces 
of the etafon. Comparing Eqs. ^ and (|T7)) . we can see 
Elf + Cib = Ci and = 62- Note that eit would simply 
become ri if nc = 1, which represents the fact that the 
light transmitting through the front surface coatings also 
probes the expansion of the coatings due to the difference 
of the refraction index of the coatings from that of the 
vacuum. 



of Eq. (|39|) . becomes different: 

[thick layer model] 
-A' 



EL, — 



A' + 2/i' 
A + 2^ 



^fcAn(0)Jo(fcr) 
XeiFopo 



Ai(3A + 2^j) /i(3A + 2^)_ 



(49) 



3.2. Elastic motion of the substrate 

Most part of the calculation process for the thick-layer 
model is common to the process for the thin-layer model 
that we have shown in Sec. 12.21 First we solve the elas- 
tic equation of the ctalon substrate and then extend the 
solution to the coatings. There are three ej coefficients 
in the thick-layer model, but we need only two weight- 
ing coefficients to be multiplied to the imaginary forces 
on the front and back surfaces of the ctalon substrate, 
and as a result the weighting coefficients are the same as 
those in the thin-layer model so that the solution of the 
elastic equation is exactly the same. It is only the expan- 
sion of the front surface coatings that makes a difference 
between the two models. 

In using the Fluctuation-dissipation theorem, the 
thick-layer model says that the imaginary force eifi^o is 
applied on the front surface coatings and the imaginary 
force eib-Fo is applied directly to the substrate. However, 
due to the thinness of the coatings (even in the thick- 
layer model), the force applied on the coatings is directly 
transferred to the substrate. The total force applied to 
the substrate is therefore eifi^o + eib^o = ei^o, which is 
exactly the same as the imaginary force on the front sur- 
face coatings in the thin-layer model. The force applied 
to the back surface is apparently €2fFo = £2-^0. 

The calculation process is common up to Eq. p7p . but 
one of the boundary conditions ([55]) becomes different: 



[thin layer model] 

nAr)=T.zir, O)--eiFop(r) 
[thick layer model] 

T^,(r) = T,,(r, 0) -t- eii,Fop{r) = -euFopir) 



which tells us that the expansion of the coatings is in 
most cases slightly overestimated in the thin-layer model. 
Accordingly, the stress tensor element Ezz, the last line 



4. RESULTS 

According to the non-zero correlation between the mo- 
tions of the front and back surfaces, the power spectrum 
of total coating thermal noise contains not only the power 
spectrum of the motion of each surface, S'i(ri) and S2{i^), 
but also the cross spectrum Sui^): 

S^in) = ejSi + 4S2 + 2eie2Si2 ■ (50) 

However, if we use Si and S2, which are respectively the 
power spectrum of noise caused by the thermal energy 
in the front surface and back surface coatings, the total 
ctalon noise spectrum Sx can be expressed without the 
cross spectrum as 

Sxi^) = (e? + ?7^e2 + 2X1 ei ^^2) Si 

+ (^'e? + C'e2 + 2X2 77eiCe2)^2 . (51) 

While the thermal motions of the surfaces can be cor- 
related, the origins of the fluctuations must be indepen- 
dent. Here the coefficient rj represents the mechanical 
transfer from one surface to the other. The coefficient 
C results from the difference of the application points of 
the probe force, and the value is slightly different be- 
tween the thin-layer model and the thick-layer models. 
With our parameters, however, ( is about 80 % in either 
model. The mechanical transfer 77 can be obtained via 
noise calculations with lossless back-surface coatings. In 
this case it equals the noise ratio, when the same force 
is applied separately onto the back and the front of the 
etalon. If we repeat this calculation with lossless front- 
surface coatings, we obtain the slightly different expres- 
sion ry/C, showing the breaking of symmetry in our model. 
The coefficients xi E^nd X2 represent the correlation of the 
motion generated from the same dissipation source but 
appearing on the different surfaces. See App.|A]for more 
details. 

The comparison of the first two terms for the back- 
surface coatings, i.e. the ryei-term and the ^e2-term, is 
shown in Fig. [3l Here we assume the proposed end mir- 
ror of the Einstein Telescope (see Table |l|. The front 
surface has 3 doublets and the back surface has 17 dou- 
blets of coating layers, which makes the reflectivity of the 
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etalon as high as that of a single mirror with 19 doublets; 
" N doublets" means that there are N quarter- wavelength 
layers of tantala, N— 1 quarter- wavelength layers of silica, 
and one half- wavelength silica cap layer. For the assumed 
mirror thickness of 30 cm, the mechanically transferred 
motion (r^ei-term) is larger than the optically transferred 
motion (Ce2-term). The measurement frequency is fixed 
at 100 Hz. 



optically transferred 
mechanically transferred 
total (back surface TN) 




— I r- 

0.2 0.4 
Thickness (m) 



0.6 



0.8 



—1 

1.0 



FIG. 3: Thermal noise originating from the back-surface coat- 
ings appears via two different paths, (i) The motion is mea- 
sured by the probe on the back surface (optically transferred), 
(ii) The motion transfers to the front surface through the 
substrate and is measured by the probe on the front surface 
(mechanically transferred) . 



TN of a conventional mirror 

TN of front surface of etalon 

TN of back surface of etalon 

— TN of etalon (total) 




— I r- 

0.2 0.4 
Thickness (m) 



0.8 



1.0 



FIG. 4: Thermal noise of etalon depends on the thickness of 
the mirror. The total noise level of the etalon can be a factor 
of ~ 1.7 times better than that of the single mirror with the 
same reflectivity when the mirror thickness is 30 cm. 



Figure |4] shows the total thermal-noise level of the anti- 
resonant etalon. Here we use the thick-layer model to 
calculate the thermal-noise level, but the difference be- 
tween the models is less than 1 %. The noise level of 
the etalon is a factor of ~ 1.7 times better than that of 
the single mirror with the same reflectivity. We could 
reduce the noise level further by u sing a thicker etalon. 
In this case, thermorefractive noise [13| will increase, but 
it turns out to be a few orders of magnitude smaller than 
coating thermal noise with our parameters |14l |. Increas- 
ing the thickness to 40 cm, we can further reduce the 



noise level by ^ 37 %. 



5. DISCUSSIONS 
5.1. Cross check with FEA 

In order to verify the results derived in the last sec- 
tion, we have performed the numerical calculation using 
a finite element analysis (FEA). With the help of the 
COMSOL program package [16| we applied two virtual 
Gaussian pressures on the front and back surface and cal- 
culated the elastic response, i.e. the strain energy distri- 
bution. In accordance with the thin layer model of Sec. [21 
firstly the elastic response of the substrate without any 
coatings was determined. Due to the small fraction of 
coating compared to the substrate, this approximation 
will be applicable. We carefully computed the error due 
to this approximation to be well below one percent. Sec- 
ondly the strain energy in the coatings results from the 
consideration of the boundary conditions between sub- 
strate and coating as is shown in Eqs. ([55)1 and (|42p . The 
numerical code was successfully tested on the analytical 
model of Brownian bulk noise given by Liu and Thorne 
[lo| . Thus, it can be regarded as an independent tool to 
confirm the semi-analytical derivation. 

We applied both results on a fused silica test mass char- 
acterized in table HI A comparison of the elastic energy 
density in the different coatings is given in table El Both 
the analytical and FEA results match within 5 %. Due 
to the Fluctuation-dissipation-theorem a coincidence of 
dissipated energy will lead to the same level of Brownian 
noise for both approaches. These results clearly confirm 
the analytical calculations. 

Additionally the FEA calculation provides a good mea- 
sure to estimate the accuracy of the quasistatic analyti- 
cal calculation. Increasing the vibrational frequency will 
lead to a failure of the quasistatic approach, typically 
starting around the lowest resonant frequencies of the 
test mass. With the help of the saiTiple's layer stack ge- 
ometry one can calculate the total dissipated energy cor- 
responding to the expected Brownian noise. Numerical 
results of the frequency behavior of the total dissipated 
energy are given in Fig. [51 They show a deviation from 
the quasistatic model of less than 10 % for frequencies 
lower than 900 Hz. 



5.2. Analytical cross check on the thin mirror limit 

Another way to verify the result in a more intuitive way 
is to compare the thermal-noise level when the mirror is 
so thin that the mechanical transfer from one surface to 
the other is almost direct (i.e. rj 1). 

The strain and stress tensors are described by the sum 
of the Besscl's functions and zcroth- and first-order poly- 
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mirror radius 


31 cm 


mirror thickness 


30 cm 


beam radius 


12 cm 


(j) of silica layers 


4 X 10"^ 


(p of tantala layers 


2 X 10"* 


number of silica layers'*^ 


2 on front, 16 on back 


number of tantala layers 


3 on front, 17 on back 


amplitude reflectivity of the front coating 


n = 0.834 


amplitude reflectivity of the back coating 


r2 = 1 - 1.36 X 10"^ 


force coefficient of the front 


ei = 0.869 


force coefficient of the back 


£2 = 0.131 


Young's modulus (silica) 


7.2 X 10^° Pa 


Young's modulus (tantala) 


1.4 X 10" Pa 


Poisson ratio (silica) 


0.17 


Poisson ratio (tantala) 


0.23 


refractive index (silica) 


1.45 


refractive index (tantala) 


2.035 


temperature 


290 K 


Boltzmann constant 


1.38 X 10"^^ J/K 


wavelength of light 


1064 nm 



TABLE I: Parameters used in this paper. A half-wavelength silica layer was placed on top of these quarter-wavelength 
doublets for protection. 



position 


material 


U [J/m] FEA 


U [J/m] analytical 


front 
front 
back 
back 


silica 
tantala 

silica 
tantala 


19.0 X 10~" 
21.7 X 10"" 
3.15 X 10"" 
6.69 X 10"" 


19.1 X 10"" 
21.9 X 10"" 
3.25 X 10"" 
6.64 X 10"" 



TABLE II: Energy density stored in the coating layers of a Khalili cavity in comparison to the analytical results at 100 Hz. 



S5 10" 



-a 10" 



1 10- 

W 
W 

b 10" 

10" 













' 2 3 4 5 6 7 8 


' 2 3 45S789' 



10' 

Frequency (Hz) 



FIG. 5: Frequency dependence of coating dissipation for 
Brownian noise. The mirror parameters are shown in Table U 



nomials of r and z. The coefRcients of the functions and 
the polynomials are given by the combination of am , , 
"fmi Co, and ci. Using an approximation h <^ a (or 



kh <fC- 1), these coefScients become 



Fpm{\ + 2ii) (ei +62) 



Ir, 
5r, 



Co 



Cl 



k^i{l+iJL) 2kh ' 
Fpm{X + 2iJ.) (ei -Hea) 
+ n) 2kh ' 

Fpm (ei + £2) 
2fc(A + n) 2kh ' 

Fpm (ei + £2) 
2fc(A + /i) 2kh ' 

Jo(fca)p, 



6a^ , , 



12a2 



ka 
Jo{ka)pr, 



(52) 
(53) 
(54) 
(55) 
(56) 

(57) 



One can see that ci and £2 appear in the coefficients in 
the form of ei -I- £2, which is actually the reflectivity of 
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the system: 



ei + £2 



ri + 7-2 
1 + rir2 



(58) 



The thermal-noise level of a conventional mirror is given 
by the same equations with ri ~ 1 and r2 = 0. The 
results of this case match with a former calculation of 



thermal coating noise of a conventional mirror [11 \ 



5.3. 



Possible noise cancellation in a resonant 
Khalili etalon 



In the case of the etalon in this paper, noise contribu- 
tions from different dissipation sources (front- or back- 
surface coatings) through different paths (optical and me- 
chanical transfer) add up to make total noise. If the light 
is not on anti-resonance in the etalon like in our case but 
on resonance in the etalon, the two noise contributions 
will not add up but partially compensate so that the total 
noise level could be lower than the anti-resonant etalon. 
The eo terms become 



s... 



- - optically transferred 

-- mectianically transferred 

— total (front-surface TN) 
anti-reso Etalon (front) 



1 1 1 1 1 

0.2 0.4 0.6 0.8 1.0 
Thickness (m) 



- - optically transferred 

-- mechanically transferred 

- total (back-surface TN) 

- ■ anti-reso etalon (back) 



0.2 0.4 
Thickness (m) 



FIG. 6: Thermal noise of an etalon with the light being reso- 
nant in the substrate. Noise contribution from the dissipation 
caused in the front-surface coatings (top) and in the back- 
surface coatings (bottom). Part of the motions probed on the 
front and the back surfaces partially compensate but the total 
noise level is higher than the anti-resonant etalon. 



reso 
^If 



ri nc(l - rl)r2 



1 - rir2 1 - rir2 
?ie(l - rl)r2 ^ Usjl - rl)r2 
(1 - rir2Y 



1 - 
ns{l 



rir2 



r2 



t2f 



(l-rir2)2 



(59) 



when the light is on resonance in the etalon. The re- 
flectivity of the etalon is equal to e^^l° + e^-^^" + e^2^° = 
(ri — r2)/(l — rir2), which is lower than ci -\- £2- With 
the parameters we used in the paper, the power trans- 
mittancc of the resonant etalon is 91 ppm while that of 
the anti- resonant etalon is 2.5 ppm. Note that the thick 
layer model is suitable in the case of the resonant etalon. 

Figure [6] shows the noise level of the resonant etalon. 
In the top panel, the origin of the fluctuation is the front- 
surface coatings, and in the bottom panel, the origin is 
the back-surface coatings. In each panel, the total noise 
level of the resonant etalon is given by the coherent sum 
of the mechanically transferred motion and the optically 
transferred motion. Since the fluctuations through dif- 
ferent paths partially compensate, the total noise level 
is lower than the noise level of the fluctuation through 
a single path. Due to the fact that the light circulates 
and probes the fluctuations more in the resonant etalon, 
however, the noise level is higher than that in the anti- 
resonant etalon in most cases. There could be a better 
parameter set to realize more significant improvement 
by the compensation, but we shall leave this for a future 
work. 



5.4. Multi layer model 

In this paper, we approximated the multi-layer coat- 
ings as a single thick layer of silica or tantala. This 
mono- layer approximation is used in Ref. Il|fl5j. Ref- 
erence [17[ explains that the stress tensor and the strain 
tensor of each layer are given by those of interface be- 
tween the substrate and the neighboring layer if we as- 
sume the stress is equal within each thin layer, thus the 
square-sum of the elastic energy in each layer is equal to 
that of a single thick layer. This approach is based on the 
assumption that the light reflects on the flrst surface of 
the coatings and the light propagates without reflections 
in the coatings. While this assumption is a good approx- 
imation, the light is actually reflected on each boundary 
between the layers. Each displacement 5xk of the fc-th 
boundary of the layers contributes to the phase shift of 
the reflected light. The total displacement 5x read out 
by the reflected beam is given as 



5x 



k 



(60) 



where is the coefficient that represents a contribution 
of each boundary. With this approach, recently demon- 
strated in Ref. [18[, the power spectrum of coating ther- 
mal noise for a conventional mirror can be corrected by a 
few percent. In the case of an etalon, the correction of the 
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power spectrum is more than the conventional mirror: 

'S'mono(^) 'S'nmlti 



0.10 



(61) 



where S',nono(^^) is calculated with the mono-layer ap- 
proximation (used in this paper) and S'nmiti(^^) is calcu- 
lated with this approach with multi-layers. 



6. SUMMARY 

In this paper we have shown a semi-analytical method 
to calculate coating thermal noise of an etalon, which 
can be used to lower the total thermal-noise level with 
the same reflectivity as a conventional mirror if we put 
a few coatings on the front surface and more on the 
back surface {Khalili etalon). We used the Fluctuation- 
dissipation theorem with probes on both sides of the 
etalon so that the thermal motion due to the same dissi- 
pation source can be summed up coherently. The result 
revealed a strong dependence of the noise contribution 
on the thickness of the etalon. In fact, the noise level of 
the etalon can be as low as that of a separate 2-mirror 
system, if the etalon is sufficiently thick. The result was 
verified with numerical calculations by a finite-element 
analysis code. The deviation was less than 5 % below 
mechanical resonance of the mass. We also checked if it 
is possible to compensate thermal motions caused by the 
same dissipation by changing the resonant condition of 
the light in the etalon. The compensation is possible but 
the noise level actually increases due to the circulation of 
the light in the etalon. 

The Khalili Etalon has clear application, for example, 
in a future gravitational-wave detector. The method we 
developed can be used also to calculate thermal noise of 
partially transmissive optics like a beamsplitter. In the 
case one uses a separate- type Khalili cavity [l| , instead of 
the etalon, our method is useful to calculate the influence 
of the back surface of the first mirror. 

There are several additional factors to be investigated 
before the Khalili Etalon is ready to be implemented in 
a future gravitational-wave detector. Several other kinds 
of thermal noise will be added by the use of the etalon. 
The etalon parameters should be optimized according to 
the magnitude of other thermal noise. From the exper- 
imental side, thermal lensing and light scattering might 
be a problem. These issues are out of scope of this paper 
but the investigation will be done in the near future [14 \ . 
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Appendix A: Mechanical transfer 



In this appendix we derive Eq. (|5ip with the parameter 
of mechanical transfer rj. We start with the application 
of forces eiFq to the front of the front coating and €2^0 
to the back surface of the substrate {Si and 62 in Fig. 
[1]). Due to the small spatial dimensions of the coatings 
we can assume stress and strain to be independent from 
the cylindrical axis in the coatings. This refers to a situ- 
ation where the forces eiFo and £2^0 £^re applied to the 
front and back surface of the substrate, respectively. In 
the approximation of thin coatings, the elastic response 
of the Khalili etalon is calculated in this paper by ne- 
glecting their influence. Then, a force acting on the front 
affects the back surface in exactly the same way as a force 
applied onto the back affects the front surface. Finally, 
the clastic energy in the coatings is obtained via tran- 
sition conditions at the substrate coating boundary (see 
Eq. 1^ and (USD). 

The linearity of the elastic equations allows us to divide 
the situation with two applied forces into two problems 
with only one force acting. This results in the following 
expressions for the energy density at position i due 
to force j. If only force ei_fo is acting on the front we 
obtain 



ii(i) = ^(eiFo)^ ,M2(i) = A{r]eiFaf 



(Al) 



The same considerations for the force on the back 62^0 
leads to 



wi(2) =A{r]e2Fof ,U2(2) ^ A{Ce2FoY 



(A2) 



The coefficient rj represents the mechanical transfer of a 
force on one side of the etalon that increases/ decreases 
the elastic energy of the coatings on the other side of 
the etalon. The symmetry of this problem is represented 
by the identical coefficient A. But note the parameter 
C breaks the symmetry of the system. This parameter 
is introduced due to the fact that the stress component 
Tzz is assumed to be zero in the back coating (i.e. the 
force €2^0 is acting on the boundary between substrate 
and coating) in contrast to the front coating. 

Examining thermal noise of the etalon demands knowl- 
edge of the dissipated energy in the coating, which de- 
pends on thickness and mechanical loss (f). We combine 
both effects into one averaged parameter 0j for each coat- 
ing i. We arrive at the following expressions for the dis- 
sipated energies AC/j, 

AC/i(i) = ^(eiFo)20i ,AC/2(i) =^(77eiFo)202 , (A3) 
At/i(2) = ^(7762^0)201 ,AC/2(2) =^(Ce2Fo)'02 • (A4) 
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To combine these two results to the case of energy loss 
with two applied forces leads to an interferometric term. 
Stress and strain behave linearly with respect to the am- 
plitude of the applied probes. But due to the spatial sep- 
aration of the probe forces the associated energies will 
not add linearly. A close examination of the problem 
reveals an analogy to the model of the classical optical 
interferometry and finally yields 

AL/i =A{e\+ rj\l + 2xir?eie2) F^'^, , (A5) 
AC/2 = A {v'el + e4 + 2x2??Ceie2) F^^^ . (A6) 

In the above expressions Xi and X2 determine the 
strength of the superposition. Due to the different con- 
ditions for Tzz in the front and back coatings the calcula- 
tion of energies differs leading to a different parameter x 
for the total energy of the front and back coating. Also 
note the relation — 1 < Xii X2 < 1 that limits the inter- 
val of possible energy values. The exact values depend 
on the solution of the elastic equations for the Khalili 
etalon. 

We want to compare this result to the situation of a 
conventional mirror. For this purpose we virtually re- 
move the back coating (02 — 0) and only consider the 
front coating of the etalon. This leads to the case of only 
a single force on the front surface (ei = 1 and €2 = 0). 
Because our calculation only considers the response of 
the substrate, we can use Eq. (|A5[) for this case, too. 
This yields the total dissipated energy 

AL/i.cM - AF^^^ . (A7) 

The same considerations hold for a conventional mirror 
only consisting of the back coating. With the help of 
Eq. (|4]) we can replace the dissipated energy by the noise 
spectrum and obtain Eq. ()5ip : 

Sxi^) = (e? + + 2xi ei ^^2) Si 

+ W4 + C4 + 2X2 ??ei (£2) S2 . 



S'l and S2 characterize the noise spectra of conventional 
mirrors possessing the etalon's front or back coating 
stack, respectively. 



From Eq. (jA5[) it is clear how to determine the co- 
efficient f]. Integrating the energy on the front coatings 
for two situations is sufficient: Calculate the dissipated 
energy of the front coating AJ/i for a force applied (a) 
on the back surface (ei = 0, £2 = 1) and (b) on the front 
surface (ei = 1, €2 = 0). Dividing the result of (a) by 
the result of (b) gives the factor rf' . If we repeat this 
procedure for the energy dissipation in the back coating 
we obtain 77^ /C^- 
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